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ABSTRACT 



We observed the Mars Trojan asteroids (5261) Eureka and (101429) 1998 VF31 
and the candidate Mars Trojan 2001 FR127 at 11.2 and 18.1 microns using 
Michelle on the Gemini North telescope. We derive diameters of 1.28, 0.78, and 
<0.52 km, respectively, with corresponding geometric visible albedos of 0.39, 
0.32, and >0.14. The albedos for Eureka and 1998 VF31 are consistent with the 
taxonomic classes and compositions (S(I)/angritic and S(VII)/achrondritic, 
respectively) and implied histories presented in a companion paper by Rivkin 
et al. Eureka's surface likely has a relatively high thermal inertia, implying 
a thin regolith that is consistent with predictions and the small size that we 
derive. 
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1 Introduction 



It has been known for nearly a century that a large population of asteroids, 
known as Trojan asteroids, exists in a 1:1 resonance with Jupiter. Neptune is 
known to have five Trojan asteroids, but the only rocky planet with known 
Trojan asteroids is Mars. At present, the number of confirmed Mars Trojan 
asteroids is four, with a handful of other candidates. All but one orbit in the 
L5 (trailing) zone, with 1 in the L4 (leading) zone. 



Mars Trojans can be dynamically sta ble over the age of the Solar System at in- 



clina tions between 12 and 40 degrees (iTabachnik and Evand . 



1999 



SchoU et al. 



20051 ): all known and candidate Mars Trojans reside in this dynamically sta- 
ble region of phase space and may therefore be primordial objects. If so, these 
bodies represent leftover planetesimals from the formation of Mars. It is there- 
fore int eresting that t here i s evidence that suggests a diverse history for these 



bodies. 



Rivkin et al 



( 120031 ) obtained visible spectra of the three largest Mars 
Trojans — (5261) Eureka, (101429) 1998 VF31, and (121514) 1999 UJ7 — and 
found that Eureka and 1998 VF31 are likely Sa- or A-class asteroids, whereas 
1999 UJ7 is probably an X-class asteroid. These differing compositions suggest 
that thes e asteroids ca.nnot h ave all formed in the same protostellar disk envi- 



ronment. 



Rivkin et al. 



(j2007l ). a companion paper, extend the previous paper, 
finding that Eureka is angritic (igneous, from an oxidized, carbonaceous chon- 
dritic precursor), whereas 1998 VF31 is likely a primitive achondrite (from a 
reduced origin). The Rivkin spectra represent the total published knowledge 
of the physical properties of Mars Trojans. Clearly, additional data are needed 
to characterize these unique objects and unravel the population's history, since 
Mars Trojans may represent the only known planetesimals that formed interior 
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to the asteroid belt, and are the closest approximations to Earth's building 
blocks currently known. 



We have obtained the first thermal infrared measurements of three Mars Tro- 
jan asteroids: (5261) Eureka, (101429) 1998 VF31, and 2001 FR127 (this last 
being unconfirmed as a Mars Trojan due to lack of long-term dynamical in- 
tegrations). Here we present our data (including some ancillary visible wave- 
length data); thermal modeling; and the resulting derived diameters and albe- 
dos. We briefly discuss the implicatio ns of our resu l ts, inc luding their general 



agreement with the taxa presented in iRivkin et al. 



mm . 



2 Recovery of 2001 FR127 



2001 FRi?7 was d iscovered in March, 2001, and recovered two weeks later 



( iTichy et al. 



2OOII ). On the basis of this arc, it was identified as a candidate 
Mars Trojan asteroid. (For simplicity, we refer this body hereafter simply as 
a Mars Trojan.) However, at the time of our Gemini observations four years 
later, the positional uncertainty for this object was around 0.67 degrees, far 
too large to be usefully targeted with Gemini/Michelle (see below), whose 
field of view is less than 1 arcminute. Before our thermal infrared observations 
could be carried out, therefore, 2001 FR127 needed to be recovered to reduce 
its positional uncertainty. 

We carried out a wide-area recovery program in May and June, 2005, using 
90prime, the prime focus camera on the University of Arizona/Steward Obser- 
vatory 90-inch (2.3-meter) Bok Telescope on Kitt Peak. 90prime has an array 
of four 4096x4096 CCDs; the sides of the array span 1.16 deg, and the total 
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2004 ). We im- 



imaged area of the sky is 1 deg^ per pointing (IWilliams et al . 
aged a series of overlapping fields along the projected location of 2001 FR127. 
We used a modified version of the Deep Eclipt i c Survey data reduction and 



moving object detection pipeline (IMillis et al 



2002 



Elliot et al. 



20051) to 



search for 200 1 FRi??, and found i t on m ultiple images. These new astromet- 



ric positions (IHergenrother et al. 



2OO5I ) allowed a substantial refinement of 



the orbital elements and a consequent reduction in positional uncertainty for 
2001 FR127. At the time of our Gemini observations of 2001 FR127 one month 
after recovery, the positional uncertainty was less than an arcsecond. 



3 Partial lightcurve photometry 

In order to detect or place limits on any visible lightcurves, we observed Eu- 
reka and 1998 VF31 in V band (predicted magnitudes from Horizons: 18.9 and 
20.6, respectively) on 2005 Sep 19 UT with the facility 2K CCD (CCD21) 
imager on the Steward Observatory Mt. Bigelow 61-inch (1.54-meter) Kuiper 
Telescope. The night was not photometric, so we cannot independently deter- 
mine the V magnitude of the asteroids. Each asteroid was easily detected in 
a series of 60 second exposures. 2001 FR127 was also attempted, but was too 
faint (predicted magnitude V=22.3 from Horizons) to be detected in 90 sec- 
ond exposures. Our time baseline was quite short - around 20 minutes for 
each of the two detected asteroids. No significant variation in flux from either 
asteroid was detected (using relative photometry to three comparably bright 
comparison stars that are nearby in the images). For the bright Eureka, our 
non-detection of flux variation places a 3a limit on any lightcurve variation of 
0.1 mag over the 20 minute observation window. This upper limit is consis- 
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tent with the iRivkin et all ( l2003l ) measurement (from a different epoch) of a 
hghtcurve amphtude at least 0.15 mag over around 6 hours. For the relatively 
faint 1998 VF31, the 3a limit on any lightcurve variation is around 1 magnitude 
over the 20 minute observation window. Both of these non-detections are use- 
ful in eliminating the possibility of fast rotations and extreme shapes, which 



would be suggested by large flux variations over these short time windows. 



4 Thermal infrared observations 



Mars Trojan asteroids Eureka, 1998 VF31, and 2001 FR127 were observed in 
queue mode on 2005 Jul 6 (UT) at the Gemini North tele scope on Mauna Ke a 



19971 ). 



using Michelle, a mid-infrared imager and spectrometer (iGlasse et al. 
In imaging mode, Michelle has a field of view of 32"x24", with pixels 0.1 arcsec 
on a side. We used the A^' (Ac = 11.2 microns, AA = 2.4 microns) and Qa 
(Ac = 18.1 microns, AA = 1.9 microns) filters (see Table [1] for the observing 
log). The telescope was tracked at asteroid (non-sidereal) rates. 

We employed the standard chop-nod strategy in which the telescope secondary 
chops (several Hertz) and the telescope nods (few times a minute) between 
nearby (8 arcsec) pointings in order to subtract out the thermal background 
from the telescope and atmosphere. The data was processed using the midir 
package of the Gemini IRAF packagJ^. The mireduce task conveniently per- 
forms all the standard tasks of reorganizing data structures and stacking the 
chop-nod images to produce final double-differenced images. In each image 
stack for this program only a single target appears, at very close to its pre- 



^ Available at http : //www . gemini . edu/sciops/data/ dataSof tware . html 
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dieted location, so there is no confusion as to whether the measured source is 
the targeted source. The exceptions are the A^' image of 2001 FR127 and the 
Qa image of 1998 VF31, in which no sources are evident at all; we derive upper 
limits to the fluxes for these two observations, as described below. In all cases, 
the images in the "off" positions are unguided and somewhat smeared, so we 
measure photometry from only the "on" (guided) positions. Each "on" source 
has one half of the total integration time. 



Three photometric calibrator stars (see ICohen et al\ (1l999l ) and the Gemini 
web pageiil) were observed using the same observing mode (Tables [Hand [2]). 
For each standard star, we measured the total flux using aperture photometry. 
The extinction we derived from these calibration measurements is consistent 
with zero. 

The Eureka and 1998 VF31 A^' observations each were made with 24 total nods. 
For the Eureka observations, the image quality clearly degrades after the first 
12 nods, so only these good quality images were used. For these data, we use 
a small centroided photometric aperture of 0.5 arcsec that is well-matched to 
the size of the image. This small aperture reduces the sky background noise 
included in the photometric aperture, but requires an aperture correction. We 
derive an aperture correction of 1.62 by measuring fluxes from the standard 
stars with both large and small apertures. The final (aperture corrected) N' 
flux for Eureka is given in Table [2l 

1998 VF31 is somewhat fainter than Eureka, and is not visible in data from 
a single nod position. We used all 24 nod positions for this target, double- 
difference combined into three intermediate images of eight consecutive nod 



http : //www . gemini . edu/ sciops/ instruments/mir/MIRStdFluxes . html 
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positions each. This allows us to make three independent measurements of the 
asteroid's brightness. The final (aperture corrected) N' fiux for 1998 VF31 is 
given in Table [21 

At Qa, Eureka is not visible in the data from individual nods nor easily visible 
in partial sums of the data. Consequently, we combine all the data into a single 
final double-differenced image. We used a small aperture of 0.7 arcsec, with 
an aperture correction of 1.23 (derived from the standard stars). The final 
(aperture corrected) Qa fiux for Eureka is given in Table [2J 

We did not detect 1998 VF31 at Qa band and did not detect 2001 FR127 at A^' 
band. (Because 2001 FR127 was not detected at A^', no measurement of that 
asteroid at Qa was attempted.) To determine upper limits on the fiuxes for 
these non-detections, we implanted scaled copies of the A^' 1998 VF31 point 
source in the 2001 FR127 A^' image and scaled copies of the Qa Eureka point 
source in the 1998 VF31 Qa image. We measure the fiuxes for the faintest 
implanted objects that are detected to set the upper limits for these two 
observations. These limits are roughly consistent with three times the sky 
noise (Table [2]), as expected. 

Eureka is the brightest asteroid source in our program, and we look to it to 
characterize the errors in our measurements. We have 12 independent mea- 
surements of Eureka at A^'. The scatter in these fiux measurements is compa- 
rable to the variation in the sky background (both around 12%). We therefore 
conclude that, for measurements where we have only few or one fiux measure- 
ment and hence cannot characterize the scatter in the measurements, we can 
estimate the error in our measurement from the sky variability. We use this 
technique for both the 1998 VF31 A^' and Eureka Qa measurements. 
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We convert Michelle counts to flux density units (mJy) by using the photo- 
metric calibrator stars (Table [2]). We combine the measurements from all the 
calibrator stars to deflne a single calibration factor for each bandpass because 
the matches in time and sky location between calibration stars and aster- 
oids are not very good. At N', the calibration factor is uncertain at the 5% 
level, as estimated from the scatter among the calibration factors from the 
three standard star measurements; at Qa, the uncertainty, derived the same 
way, is around 10%. We add these errors in quadrature with the errors in our 
measurements to derive the flnal errors reported in Table [21 

Our targets have temperatures around 250 K (as calculated in our thermal 
modeling, described below). Thus, color corrections can be important, as these 
asteroids are substantially colder than the calibrator stars (~4000 K). The 
isophotal wavelengths (wavelength at which the flux density from the object's 
spectrum equals the average flux density calculated by integrating over the 
fllter proflle) for Ceres (217 K) in Michelle N' and Qa are 11.52 microns and 
18.26 microns, respectively (K. Volk, priv. comm.). We adopt these Ceres 
isophotal wavelengths for our observed asteroids as our derived temperatures 
are quite close to the Ceres temperature used. We use these isophotal wave- 
lengths in our thermal models described below, although these shifts from 
nominal wavelengths are relatively small and, compared to the errors in our 
measurements, unimportant. 



5 Thermal modeling and model results 



We interpret the observed thermal emission from our tar g ets us ing the Stan 



dard Thermal Model (STM; e.g. 



Lebofskv and Spenceil (119891 )). The STM 



10 



assumes a non-rotating (or zero thermal-inertia) spherical asteroid: dayside 
temper atures are in equilibriu m with sunlight, while the nightside temperature 



IS zero. 



Lebofsky et al. 



(119861 ) found that the thermal emission from asteroids 
frequently has a higher color temperature than would nominally be predicted 
under the STM assumptions given above, and introduced an empirical param- 
eter, ?7, that allowed them to simultaneously model the elevated color temper- 
ature of the thermal emission and the (known) size of their targets (Ceres and 
Pallas) . The canonical value for this beaming pararn e ter r? is 0.756, but recen t 



studies {e.g., 



Harris ( 



Stansberry et al 



19981) 



Delbo and HarrisI (120021 ) 



Fernandez et al. 



mm 



(120071 )) find that r] can be significantly larger. 



We rely on JPL's Horizons ephemeris service for distance, phase angle, and 
absolute visual magnitude {Hy) information (Table [H Table [2]). We relate 
diameter, albedo, and Hy through 



D 



1329 



X 



where vy is t he visible geometric albedo and D is the diameter in kilometers 



( iHarris 



19981 ). We use a thermal phase coefficient 0.01 mag/deg. We assume 
standard scattering behavior for the surface in the visible, resulting in a phase 
integral q = 0.39. 

Our data on Eureka allow us to determine the albedo pv, diameter D, and 
beaming parameter rj because Eureka is strongly detected at both A^' and 
Qa (in combination with Hy, this gives three measurements and three un- 
knowns). Models with t] = 1.3 give the best fit to the observed A^' — Qa color 
(Figured]), but t] values in the range 0.57-2.45 are all formally consistent with 
the data and error bars given in Table El Using this range of rj, we derive a 
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diameter for Eureka of 1.28lo;3| km. The corresponding alb edo, when 



the di- 



ameter uncertainties and aii uncertainty of 0.3 mag for Hy ( iJuric et al. 



Romanishin &: Teglerl . 



2002 



20051 ) are included, is 0.39lo;23, a range of acceptable 
albedos that is so large as to be nearly useless. We take two approaches to 
generating more useful error bars. 



Our first approach relies on results from other Solar System observing pro- 
grams. Modeling of various Solar System observations has sh own that the 



beaming parameter r? is unlikely 


2003; 


Delbo et al.. 


2003, 


2007: 



y to have a value larger t han ~1.8 (IFernandez et al 



Stansberry et al. 



20071 ). Requiring r] < 1.8 re- 



duces the upper bound on Eureka's size and consequently the lower bound on 
the albedo, giving a diameter of 1.28lo;34 km and an albedo of 0.39to;ig (in- 
cluding the unce rtainty in Hy). Furth e rmore, Eureka's albedo is unlikely to be 



large r than 0.5 (IHarris and Lagerros 



2002 



Delbo et al. 



2003 



Rivkin et al. 



20071 ). which constrains the upper bound of the albedo and the lower bound 



of the diameter, giving 1.28lo;29 km and and an albedo of 0.39lo; 



0.11 
18- 



Our second approach simply assigns t] = 1.3 (the best fit for r]) rather than 
using a range of r] values. This gives a diameter of 1.28io;og km and a resulting 
albedo of 0.391qj3 (with the uncertainty in Hy still included). We take these 
latter results as our derived best fits (Table [2]), but emphasize that diameter 
and albedo both depend on choice of 1], whose uncertainty is not captured in 
the error bars for this best fit, and that alternate assumptions produce different 
"best" results (as in our first approach). However, we note that in all cases the 
derived (and best-fit) albedo is consistent with the asteroid's S (I) taxonomic 



class and interpreted angritic composition (IRivkin et al. 



20071 1. 



Drawing upon the success of the STM in fitting the Eureka data, we apply 
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the STM with r] = 1.3 to the data for 1998 VF31 and 2001 FR127; the re- 
sults are shown in Figure [T] and Table O This extrapolation of the STM to 
these other two bodies is warranted both because the sizes and temperatures 
of these three asteroids do not differ significantly, and because the available 
data for 1998 VF31 and 2001 FR127 are so sparse that they do not justify a 
different model. Because we have only a single mid-infrared data point for 
1998 VF31, we solve for the albedo and diameter (Table [2]) simply by finding 
the STM solution that passes through the A^' measurement. The errors on di- 
ameter give the range of solutions that are consistent with the la photometric 
error bars; the albedo uncertainties include both the diameter uncertainty and 
0.3 magnitudes uncertainty in Hy- For 2001 FR127, we solve for the maximum 
radius and minimum albedo by finding the STM solution that passes through 
the N' upper limit measurement. The lower limit on albedo here includes an 
uncertainty of 0.3 magnitudes for Hy- 



5.1 Eureka's thin regolith 



Spencer et al\ (119891 ) define the dimensionless thermal parameter 6, which 
gives the ratio of the characteristic radiation timescale to the diurnal (rotation) 
timescale: 



e 



ea [(1 - A) S^ 



where F is the thermal inertia of the asteroid (units of J/m^/K/sec^/^); to is 
the rotational frequency {u = 2*tt/P, where P is the rotational period); R is 
the heliocentric distance to the asteroid, in AU (Table [1]); e is the emissivity 
(we use e = 0.9); a is the Stefan-Boltzmann constant; A is the bolometric 
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albedo, which is j9y x g = 0.137; and 5*1 is the solar constant at 1 AU (we use 
1366 W/m^). For small asteroids with rj = 1.3 and phase angle ~30 degrees, 



r is expected to be 200- 400 J/mVK/sec^/^ ^nd 6 is likely to be 1-2 JPelbo 



2004 



Delbo et al. 



2003) 



Rivkin et al\ (120031 ) show a partial lightcurve for Eureka that suggests that 
Eureka's rotation period may be ~10 hours. This period gives G of 1.6-3.2 
(for the above range of F), in good agreement with the expected value. Thus, 
Eurek a's thermal p roperties appear to be consistent with theoretical expecta- 



tions (iDelbo 



objects (NEOs) flDelbo et al. 



20041) as well a s with the results for similarly-sized near Earth 



20071 ) 



In the STM, dayside temperature goes as ri~^/'^^ so r/ < 1 results in warmer 
emission, while rj > 1 results in cooler emission. Relatively cooler dayside 
temperatures will be produced by a high thermal inertia surface since some of 
the thermal emission will occur on the nightside. Our result that rj is greater 
than unity is therefore consistent with our argument that F (thermal inertia) 
is relatively large. Both conclusions are co nsistent with a relative dearth o f 



regolith, as predicted for bodies of this size ( 



and with the small size we derive {e.g., 



Binzel et al 



Delbo et al. 



2004 



Cheng 



20051) 



fl2007h l We note that rel- 



atively large rj could also indicate a relatively smooth surface (on macroscopic 
scales), which again could imply a relatively thin regolith. 



6 Summary and discussion 

We observed the Mars Trojan asteroids Eureka, 1998 VF31, and 2001 FR127 
at N' (11.2 microns) and Qa (18.1 microns). Using the Standard Thermal 
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Model, we derive diameters (albedos) of 1.28 km (0.39), 0.78 km (0.32), and 
<0.52 km (>0.14), respectively. From several lines of argument, we conclude 
that Eureka's regolith is likely relatively thin, as predicted for a body this size. 
Eureka's thermal inertia is likely similar to that for comparably-sized NEOs. 



Rivkin et al\ (120071 ) find that Eureka is angritic and that 1998 VF31 appears to 
be an S(VII) (e.g., primitive achondrite) asteroid. In both cases, the albedos 
we derive are consistent with those taxa and implied compositions. 



2001 FR127 has a diameter (<520 m) comp arable to ve r y sm all NEOs that 



have been observe d in the mid- i nfrare d {e.g. 



Delbo et al. 



experiments {e.g., 



Ostro et al. 



( I2OO3I )) and in radar 



(120021 )). These are among the smallest objects 
studied in the Solar System. By this virtue, further physical studies (spin, 
shape, lightcurve, etc.) are interesting on their own. Additionally, comparing 
physical properties between the dynamically old Mars Trojan asteroids studied 
here and the comparably sized but dynamically young NEOs may be useful 



in understanding the evolution of these smallest asteroids. 
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Observing log 



Target 


Filter 


Obs. time 


Exp. time 


Airmass 


R 


A 


phase 






(UT) 


(sec) 




(AU) 


(AU) 


(deg) 


Eureka 


N' 


9.1675000 


60 


1.33 


1.49 


0.67 


34.4 


Eureka 


Qa 


9.5359725 


150 


1.44 


1.49 


0.67 


34.4 


1998 VF31 


N' 


10.783750 


120 


1.38 


1.54 


0.81 


37.2 


1998 VF31 


Qa 


11.152500 


150 


1.50 


1.54 


0.81 


37.2 


2001 FR127 


N' 


11.989445 


120 


1.21 


1.69 


0.93 


31.4 


HD 141992 


N' 


8.7205560 


15 


1.07 








HD 141992 


Qa 


8.8315280 


45 


1.07 








HD 156283 


N' 


11.471250 


15 


1.31 








HD 158899 


N' 


12.296390 


15 


1.46 








HD 158899 


Qa 


12.407080 


45 


1.48 









Table 1 



All observations were made on 2005 Jul 06 (UT), with the exact midtimes indicated. 
The target asteroids and calibrator stars are given in the top and bottom groups, 
respectively. The exposure times used to measure object flux are listed here. The 
total open-shutter times, including both "on" and "off" (unguided) images, were 
240 and 300 sec for asteroids and 30 and 90 sec for calibration stars at A^' and Qa, 
respectively. However, for all targets, only the "on" (guided) images were used, and 
for Eureka, the second half of the N' data was poor and not used. 
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Photometry and modeling results 



Target 


H 


V 


N' 


Qa 


Albedo 


Diameter 




(mag) 


(mag) 


(mJy) 


(mJy) 




(km) 


Eureka 


16.1 


17.51 


21.7 (2.8) 


30.8 (4.8) 


0.39t°:li 


1 90+O.O8 
-■■•^"-o.oe 


1998 VF31 


17.4 


19.38 


5.1 (0.9) 


<15.4 


0.32t°:J? 


0.78 ± 0.06 


2001 FR127 


18.9 


21.21 


<1.53 




>0.14 


<0.52 



Table 2 

H magnitudes; V magnitudes at time of Michelle observations, from Horizons; 
measured fluxes at the isophotal wavelengths of 11.52 microns and 18.26 microns 
for N' and Qa, respectively; and derived physical properties (using the modified 
STM). The uncertainty in H (and therefore V) is taken to be 0.3 mag. The la 
errors for our derived albedos and diameters are given. No observation was made of 
2001 FR127 at Qa- These STM model resuhs use rj = 1.3. We used the following N' 
fluxes for photometric calibration: HD 141992, 13.586 Jy; HD 156283, 35.096 Jy; 
HD 158899, 11.754 Jy. We used the following Qa fluxes for photometric calibration: 
HD 141992, 5.362 Jy; HD 158899, 4.650 Jy (Photometric calibration fiuxes from 



Cohen et al. 



(jl999l ) and Gemini web pages.) 
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Wavelength (microns) 



Fig. 1. Observed fluxes (plotted at isophotal wavelengths) and spectral energy 
distributions (solid curves) for the three observed Mars Trojans (as labeled), with 
STM solutions for a range of albedos and rj = 1.3. Our detections are indicated 
by plotted error bars; upper limits are shown as downward pointing triangles. The 
isophotal wavelengths for these asteroids (temperatures around 250 K) are 11.52 and 
18.26 microns for N' and Qa, respectively (K. Volk, priv. comm.). Data in the visi- 
ble (circles) are taken from JPL's Horizons service. (Error bars on the visible fluxes, 
corresponding to uncertainties in H of 0.3 magnitudes, are shown, but are generally 
smaller than the symbol size.) Based on our lack of detection of photometric vari- 
ation (lightcurve) , we assign no extra scatter to the reflected light (visible) data. 
For each case, we derive diameters from thermal fluxes (assuming rj = 1.3) and 
derive albedo from the derived diameter and known Hy- Here we show models that 
correspond to six albedos: 0.05, 0.10, 0.20, 0.30, 0.40, and 0.50, from top to bottom, 
respectively. These models assume the nominal H values for each asteroid; larger 
H values imply smaller albedos and smaller H values give larger albedos. 
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